"Cosmological" quasiparticle production in harmonically trapped superfluid gases 



O 

o 

(N 



Petr O. Fedichev^^'^ and Uwe R. Fischer^ 
^ Leopold- Franzens-Universitat Innsbruck, Institut fur Theoretische Physik, Technikerstrasse 25, A-6020 Innsbruck, Austria 
^Russian Research Center Kurchatov Institute, Kurchatov Square, 123182 Moscow, Russia 

We show that a variety of cosmologically motivated effective quasiparticle space-times can be 
produced in harmonically trapped superfluid Bose and Fermi gases. We study the analogue of 
cosmological particle production in these effective space-times, induced by trapping potentials and 
coupling constants possessing an arbitrary time dependence. The WKB probabilities for phonon 
creation from the superfluid vacuum are calculated, and an experimental procedure to detect quasi- 
particle production by measuring density-density correlation functions is proposed. 
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I. INTRODUCTION 

In a gravitational field with explicit time dependence 
in the metric, particles and antiparticles can be simulta- 
neously created by quantum fluctuations from the vac- 
uum. By the uncertainty principle, the time scale of the 
system's evolution dictates the typical energy of the par- 
ticles produced 0. The process of cosmological particle 
production, whose condensed matter analogue we shall 
consider here, is potentially relevant in the expanding 
early universe, in which phonons experience an acoustic 
geometry; as a consequence, the expansion of the uni- 
verse could generate density waves growing into galaxies 


Attention has of late focused on condensed matter 
analogs of the curved space-times familiar from grav- 
ity, primarily due to their conceptual simplicity and re- 
alizability in the laboratory |llllllllllalilI3- 

Condensed matter systems lend themselves for an ex- 
ploration of kinematical properties of curved space-times 
and, in particular, provide a testbed to study the ef- 
fects of a well-defined and controlled "trans-Planckian" 
physics, i.e. atomic many-body physics on a microscopic 
scale, on low-energy quantum effects like Hawking radi- 
ation nH and cosmological particle production. In the 
present paper, we investigate quantum fields propagating 
on effective curved space-times backgrounds, for the case 
of harmonically trapped, dilute superfluid gases with pos- 
sibly time varying particle interactions. For a perfect, ir- 
rotational liquid, described by Euler and continuity equa- 
tions, it was recognized by Unruh that the action of 
fluctuations of the velocity potential around a spa- 
tially inhomogeneous and time dependent background, 
can be identified with the action of a minimally coupled 
scalar field according to 



the summation convention over equal indices, unless indi- 
cated otherwise. The quantities g^^ are the contravariant 
components of the effective metric tensor related to its 
covariant components by g^'^g^a — ^^a, and g = detg 
is the determinant of the metric tensor. The action 
leads to the "relativistic" scalar wave equation 



1 



(2) 



In general, the effects of quantum fluctuations described 
by the quantum version of Eq. are very small and 
can hardly be observed because of finite temperature and 
dissipation effects. Therefore atomic superfiuids, where 
both extremely small temperatures and dissipationless 
fiows are possible, attract growing interest for an emerg- 
ing research field of "experimental cosmology." 

In the following, we study how various curved space- 
times can be implemented in harmonically trapped su- 
perfluid Bose and Fermi gases. As a concrete example, 
we show how de Sitter and Friedmann-Robertson- Walker 
(FRW) universes can be "re-created" in superfluid gases. 
We analyze the quasiparticle production probabilities, 
leading to a thermal spectrum in the WKB approxima- 
tion, and discuss an experimental procedure to observe 
and characterize the excitations produced. 



II. QUASIPARTICLE METRIC TENSORS IN 
HARMONICALLY TRAPPED SUPERFLUIDS 

A. Superfluid action 

The hydrodynamic, i.e. long-wavelength action of a 
trapped superfluid is generally given by 



S = / dtSx— 
2k 



d 

ot ' 



dtd^x^g'^'^dM^ . 



Here, v is the background velocity, 1/k the compress- 
ibility, and c the speed of sound of the liquid. We use 



S = 



dtd^3 



d 



P, 



P^0+|(V</))'+e(p)+l^trapP 



, (3) 



where the external harmonic potential Vfcap{x,t) = 
^ ' \{uj1{t)x'^ -I- U!y{t)y'^ + Lo1{t)z'^) is characterized by the 
three frequencies lOi, i — x,y,z. The trapping frequen- 
cies are assumed to be time dependent in an arbitrary 
manner (we can also conceive of making Lof (t) effectively 
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negative by "turning over" the potential, see section llTT|l . 
In the above action, the quantity p plays the role of a re- 
sponse or stiffness coefficient to gradients of 4>, and equals 
the total fluid density at absolute zero; the equation of 
state of the superfluid is given by the energy density func- 
tional e — e{p). (Note that we leave out an overall minus 
sign in the definition of the action S.) We generally set 
h = m = 1. The action entails the existence of a veloc- 
ity potential (j), such that the vorticity is zero except on 
singular lines, and ensures the validity of Euler and con- 
tinuity equations for the superfluid velocity v = V0. The 
identification of (j) with the phase of a complex "order pa- 
rameter" (i.e., the direction of a unit vector in the plane 
of some abstract space) leads to the quantization of cir- 
culation, because (f) is then defined only modulo multiples 
of 2tt. Finally, the above action implies the conjugateness 
of phase and density quantum variables jl5j | : 



[p{x) , (j){x')] — iS{x — x') . 



(4) 



These properties, taken together, constitute the canoni- 
cal definition of a superfluid at T = 0| . There- 
fore, Eq. ^ represents the universal action of a simple 
scalar superfluid at absolute zero made up of elementary 
bosonic or fermionic atoms, independent of a particular 
microscopic model. 

The simplest example for the equation of state is that 
of a weakly interacting Bose gas, with cb = ■^9P^, where 
g is the coupling constant, g = 47ras, with as the s-wave 
scattering length, characterizing pair collisions of atoms. 
The scattering length can be tuned using external mag- 
netic fields |13| . Another example is a two-component 
Fermi gas with attractive interactions between atoms of 
different hyperfine species Ti*] . The ground state of such 
a gas is superfluid (in the simplest version, it is the BCS 
state of a scalar superfluid with s-wave pairing), and since 
the interactions are weak, the BCS gap is small and the 
equation of state (to exponential accuracy) coincides with 
that of a free Fermi gas: ep = [{Sir'^y/^ /lO] p^/^. To con- 
sider all possible cases which have a power law density- 
dependence of the equation of state, in a generic way, we 
write 



e(p) = Ag^p-', 



(5) 



where A is a numerical constant. That is, /3 = f , 7 = 2 for 
a dilute Bose gas, and /3 = 0,7 = 5/3 for noninteracting 
two-component fermions. 

In our present context, an important quantity char- 
acterizing a superfluid is the so-called "Planckian" en- 
ergy scale, Epi, i.e. the frequency beyond which the 
spectrum of the excitations above the superfluid ground 
state ceases to be phononic and (pseudo-)Lorentz in- 
variance is broken. For a weakly interacting Bose gas 
Epi ^ gp, of order the mean interparticle interaction. 
In a BCS superfluid, Epi is determined by the BCS gap: 
Epi ^ p^/^ exp[— l/(|g|p^/^)]. The complete analogy with 
a quantum field theory on a fixed curved space-time back- 
ground given by ^ only exists if all timescales to, de- 



scribing the evolution of the superfluid, are much larger 
than the "Planck time": to '3> l/Epi. 



B. Scaling transformation for Bose and Fermi 
superfluids 

Our approach in the following is based on the so-called 
scaling transformation [l^ 120, I^J to describe the ex- 
pansion and contraction of the gas under time dependent 
variations of the trapping frequencies. It is by now a 
well-established fact that the hydrodynamic solution for 
density and velocity of motion for such a system may be 
obtained from a given initial solution by a scaling pro- 
cedure both in the bosonic j3 Qi 113 as well as in the 
fermionic case :21,]. Defining the scaled coordinate vec- 
tor ajft EiXi/bi, density and velocity are given by the 
scaling transformations p!a |: 



p{x,t) 



V 



{x,t) ^ -^xf + 4>{Xb,t) . 



(6) 
(7) 



The (dimensionless) scaling volume V = Yli in den- 
sity jnj is dictated by particle conservation. 
Introducing a new "scaling time" variable by 



dTs ^ {glg{0)f 
dt V-i ' 

we rewrite the action Q in the form 



(8) 



S = 



dTgd^Xh 



OTs 2 
+ l{p) + V"trap(a;fc,0)p] , 



(9) 



where the Ts dependent scaling factors are 



V7- 



-— (10) 
6? dr. ' ^ > 



and e(p) — e{p)\g^g(Q)] "S/u = d/dxu- The rescaled den- 
sity p has no explicit dependence (whereas it has ex- 
plicit dependence on the lab time t), and coincides with 
the equilibrium condensate density profile in the scal- 
ing coordinate Xb. Where any confusion might arise, we 
will generally designate scaling variables with a tilde to 
clearly distinguish them from lab frame variables. 

For the relation H1U|I between the scaling factors and 
hi, g to hold true, we must impose the following equations 
of motion for the scaling parameters bf. 



h+uj1{t)bi = 



(g/.g(o))^^§ 



(11) 



They need to be solved with the initial conditions 6j = 1 
and bi = 0. Note that here no summation convention is 
used in the second term on the LHS. For a sufficiently 
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large cloud, the stationary background solution can be 
found from the Thomas-Fermi density profile. It is given 
by using that /i = —d(j)/dTs equals the initial chemical 
potential, and 



di _ 
dp ^ 



Vtrap(a;fc,0) 



(12) 



The part of the action quadratic in the fluctuations is 
obtained to be 



5(2) = / dTsd^Xb 



-nSp 



(13) 



where Sp = p — po and S(f> = (f> + pr^ . The rescaled bulk 
compressional modulus (inverse compressibility) k — 
d^e/d^p does not depend on the time r^, and is iden- 
tical to g(0) in the bosonic case. After integrating out 
the density fluctuations, we obtain the effective action 
for the rescaled phase variable 



1 



2k 



- 1 1'* 



where the squared scaling speed of sound — kp{xb). 
Using the identification with a minimally coupled scalar 
field, analogous to the one performed in the second line 
of Eq. ^ , the line element in the scaling variables reads 



(15) 



The line element takes a particularly simple form for an 
isotropic superfiuid Fermi gas, where 7 = 5/3 and thus 
all Fi = 1, leading to drg/dt = b^^. We note that even for 
this simple case, the metric defined by (|15|) is not trivial, 
since Ts and t are different, and both c and k depend on 
the radial scaling coordinate rf,. We will see below that in 
the case that = 1 the scaling transformation is exact, 
and that therefore no quasiparticle creation occurs in the 
scaling variable basis, i.e. there is no mixing of negative 
and positive frequency parts in the time Tg (quasiparticle 
creation can take place in the lab frame with time t, 
though, and a lab detector will still see that quasiparticles 
are "created"). 

Identifying Scf) with $, and going back from scaling 
coordinates to laboratory frame variables, we recover the 
action with v — {bi/bi)riei, and the line element, 
which is of Painleve-GuUstrand type, reads [3 11^ 



ds^ 



[~{c^ - v'^)dt^ - 2vidxidt + dx}] 



(16) 



where = up is the squared instantaneous speed of 
sound. We now assume that space is spherically sym- 
metric, i.e. that v has a radial component Vr only, 
that Vr/c = f{r) holds, and furthermore c = c{t) is 
a function of time only. We first apply the transfor- 
mation codt = c(t)dt, where Cq is some constant (ini- 
tial) sound speed, connecting the laboratory time t to 



the time variable t. 



This results in the line element 
P)dP - 2fcodtdz + dr^ + r^dVl?. 
We then employ a second transformation c^dr = codt + 
fdr/{\ — P) [23, to bring the metric into the form 



= (c/k)[-c2(1 



ds' 



1 



fir)) 



1 



1-P{r) 



dr^ + r^d^' 



(17) 

The metric in the above form facilitates comparison with 
metric tensors in spherically symmetric space-times writ- 
ten in their standard form. E.g., if p{r) — 2M/r is 
chosen, this line element is conformally equivalent to the 
Schwarzschild metric, the asymptotically flat vacuum so- 
lution of the Einstein equations around a spherically sym- 
metric body with total mass M p^ . 



III. CREATING DE SITTER AND 
FRIEDMANN-ROBERTSON- WALKER 
UNIVERSES 

The equation of state for Bose superfluids contains the 
interatomic interaction. Therefore, by varying this inter- 
action, possibly together with the trapping frequencies, 
expanding clouds of Bose atoms allow for the simulation 
of a large set of cosmological space-times. We begin by 
discussing the so-called de Sitter universe, which is a so- 
lution of the vacuum Einstein equations characterized by 
the line element El 



ds^ 



1 



A 



dP 



1 ^ 2 



dP + pdif 



(18) 

with A being the cosmological constant = energy density 
of the vacuum. Up to the conformal factor c/g, the metric 
(fTTjl coincides with the de Sitter metric (fT^ . provided we 
require that v^/c^ = Ar^/3 and that the speed of sound 
is a constant in space and time. The speed of sound in 
the center of the cloud is time independent if: 



(? = const. 



5(i)/g(0) = b\t) . (19) 



Close to the center of the condensate, c is, in addition, 
practically spatially independent. Using Eq. Q, we find 
that A = const, provided b oc exp[At], with A = c^J A/i. 
This exponential expansion of the cloud can (asymptot- 
ically) be achieved if we turn over the potential, mak- 
ing it expel the particles rather than trapping them: 



w2(t 



00 



-A . The de Sitter horizon, where Vr 



c, 



is stationary and situated at ru — c(0)/A, which is well 
inside the expanding cloud provided A 3> wq, where ojq 
is the initial trap frequency. 

The experimental sequence leading to "condensate in- 
fiation" is schematically depicted in FigQ] We assume 
that the experiment can be done with one trapped (low- 
field seeking) and one untrapped (high-field seeking) hy- 
perfine component of the same atomic species. We start 
from a sufficiently large Bose-Einstein-condensed cloud 
at small (effectively zero) temperature with all atoms 
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being in the trapped state. Then, we transfer all the 
atoms to the untrapped state, by flipping the sign of the 
trapping potential. At the same time, we ramp up the 
interaction strength, according to condition (|19|l . using 
a suitable Feshbach resonance 13]. As a result of the 
simultaneous action of the inverted parabolic potential 
and the increasing interaction energy, the gas cloud ex- 
periences a rapid exponential expansion, representing the 
analogue of cosmological inflation. 

In fact, Eq. I|19|l defines a broad class of Bose superfluid 
effective space-times. In the present case of isotropic ex- 
pansion with bi = 5, we have t = Ts and, in scaling vari- 
ables, we obtain up to a conformal factor a Friedmann- 
Robertson- Walker metric: 



c 



(20) 



According to the above form of the metric, the quan- 
tity b plays the role of the scaling parameter not only 
in our condensed matter context, but can be interpreted 
equally well as the scale factor in the expansion of the uni- 
verse, with H = b/b the Hubble parameter. As demon- 
strated above, exponential growth of 6, with constant H, 
corresponds to exponential inflation [221 ■ The present 
setup also allows for the simulation of power law infla- 
tion milj, with b{t) cx t^. The "Hubble parameter" 
H changes for all exponents 7 inversely proportional to 
time t, H cx and the exponent 5 = 1/2 corresponds 
to a "radiation dominated" universe, while the exponent 
(5 = 2/3 corresponds to a "matter dominated" universe. 
An isotropically trapped expanding superfluid gas thus 
models an isotropic expanding universe. Generically, we 
can model anisotropic universes with H15|l . with scaling 
factors which are different in different spatial directions. 
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rf 






F = - 



lower-dimensional systems; e.g., a 1-l-lD analogue of a de 
Sitter universe can be achieved for quasi-lD excitations 
in a linearly expanding elongated Bose-condensate, with- 
out changing the interaction ^2^. Another possibility is 
to use superfluid Fermi gases. We have, in the isotropic 
case, F — 1, drg = dt/b'^, and the metric may be written 
in the form 



''dt^ + p^dn^ + p^n^dn^ 



(21) 



where we defined a scale factor /3 = 6^. Performing ex- 
periments in superfluid Fermi gases has the advantage 
that three-body losses are strongly suppressed by Fermi 
statistics ISOl. 



IV. COSMOLOGICAL PARTICLE 
PRODUCTION ANALOGUE 

Now we turn to describe the evolution of quantum fluc- 
tuations, on top of the classical (mean fleld) hydrody- 
namic solutions described above. The equation of motion 
for the phase fluctuations can be obtained after variation 
of the action (|14(l . 



K^M ( —Fi{Ts)\'bti 
K 



= 0. 



(22) 



Phrased in curved space-time language, the above equa- 
tion is the minimally coupled massless scalar wave equa- 
tion for (5/5, analogous to 0, with the metric H15|l . 
Consider for simplicity the isotropic case. 



537- 



(5/5(0))'^" 

The solution for the full quantum fleld reads 



(23) 



2Ven 



4>niXb) [anXn + alXn] > (24) 



where V is the the initial Thomas-Fermi volume of the 
cloud, the operators a„ (o|i) annihilate (create) phonon 
excitations in the initial vacuum state, and the mode 
functions Xn satisfy 



Xr. 



■Fir,)e^Xn 



0. 



(25) 



FIG. 1: Exponential expansion of a two-component Bose- 
Einstein condensate with effective spin F — 1/2, to create a de 
Sitter quasiparticle universe close to its center. The trapping 
potential is inverted by applying a radio frequency (rf) pulse, 
which transfers the atoms to their untrapped hyperfine state. 

While this experiment is feasible in principle, increas- 
ing the interaction dramatically increases three-body 
losses as well, whose total rate scales like g'^p^- This 
complication can be avoided, by switching to effectively 



The initial conditions XniTs — *■ — oo) = exp[— ie„rs] are 
selected such that Eg. (|24|l . at t — > — oo, describes the 
phase fluctuations in a static trapped superfluid in its 
ground state. In quantum fleld theory (QFT) language 
this ensures that a laboratory frame detector does not 
detect quasiparticles at t — > — oo. Hereafter we define 
the "scaling vacuum" to be the quantum state annihi- 
lated by the operators a„ , where our choice of the initial 
conditions gurantees that the initial superfluid vacuum 
and the scaling vacuum coincide a,t t — oo. 
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The case when all = 1 is remarkably special: In 
this case Eq. (1^ does not depend on the superfluid evo- 
lution, and thus the quantum state of the excitations 
remains unchanged. As we have seen above, this indeed 
happens in the case of an isotropic Fermi superfluid. An- 
other example is a 2D isotropic dilute Bose gas with con- 
stant particle interaction 19]. In these cases the scaling 
transformation is exact, both for the condensate and the 
excitations. In the language of QFT this amounts to 
the fact that there is no particle production in the scal- 
ing basis, since the scaling solution is constructed from 
eigenf unctions of the (exactly conserved) scaling trans- 
formation operator Fi (in other words the scaling vac- 
uum is protected by an exact scaling invariance which 
forbids frequency mixing). The fact that no excitations 
are produced in the scaling basis does not mean that a 
lab detector does not detect quasiparticles. The trans- 
formation from the laboratory time t to the scaling time 
Ts time is nontrivial and thus the phase of the functions 
Xn is a complicated function of the laboratory time t. In 
other words, the phase field given by Eg. 12411 . if coupled 
to a detector of the type considered in [22| , gives a non- 
zero response. This indeterminacy of the vacuum state 
finds its coun terp art in the Unruh-Davies effect in flat 
space-time 0, 0| and its curved S£ace-time generaliza- 



tion, the Gibbons-Hawking effect 

The description of a quantum field state in terms of 
particles and antiparticles is based upon the separation 
of positive and negative frequency parts. As we con- 
fine ourselves to a measurement involving the laboratory 
time variable t, this distinction is only possible if the 
asymptotic phase of x„-functions is large and sufficiently 
quickly increases as a function of t. Using a WKB ap- 
proximation to the solutions of Eq. 1)25(1 , we find 



/* dr 
^/F—ldt ^ 0. 
dt 



(26) 



The latter condition can be also called a "Trans- 
Planckian safety condition" (TP condition), since if full- 
filled it implies that an experiment in a lab frame prob- 
ing an energy scale Eq <C Epi does not require infor- 
mation about solutions of Eq. (|25|l with e„ > Ep\. For 
isotropic expansion of a 3D Bose gas, Eq. (|26|l is equiva- 
lent to divergence of J dtg^^^ /b^^"^ and is quite restrictive: 
For the FRW analogy discussed above avoiding the diver- 
gence implies, according to (|19f) . that b should not grow 
faster than linearly. The TP condition (|26|l is based upon 
the WKB approximation condition for Eq. H25|) . leading 
to the requirement F(ts) ^ (cnTs)"^ for large (here 
> means "grows faster than"). Substituting the latter 
condition into Eq. (|26|l we find that the marginal WKB 
case corresponds to a logarithmically divergent integral 
in Eq. H26|l. Thus the marginal TP case corresponds to 
the marginal WKB case and vice versa. 

The equation (|25|l is formally equivalent to scattering 
of a non-relativistic particle with energy e„ by a potential 
e^(l — F{ts)). The initial conditions correspond to a 
single particle per unit time incident on the potential 



barrier. Time dependence of the scaling factors leads to 
scattering of the particles from the incoming wave and 
at Ts —> oo the WKB solution reads: 



1 



where an is the transmission and Pn the forward scat- 
tering amplitude. The coefficients an and are related 
via the particle flux conservation condition: 



|anp-|/3„P = l. 



(28) 



In QFT language the latter condition is the Bogoliubov 
transformation normalization condition for a bosonic 
field. The number of particles detected by a scaling time 
detector at rest is measured by the absolute square |/3„p 
(which is proportional to the probability that the detec- 
tor absorbs a quantum), and therefore Nn = |/3nP can be 
interpreted as the number of scaling basis quasiparticles 
created. 

In the WKB approximation the amplitudes are con- 
nected in a simple way: 



exp[ 



j2To]ar. 



(29) 



where the inverse temperature is given by the integral 



1 



FdT, 



(30) 



and C is the contour in the complex Tg-plane enclosing 
the closest to the real axis singular point of the function 
F{ts) 0. Together with Eq. this gives 



Nn 



exp[e„/To] - 1 ' 



(31) 



i.e. adiabatic evolution of trapped gases leads to "cos- 
mological" quasiparticle creation with thermal occupa- 
tion numbers in the scaling basis. The temperature To 
depends on the details of the scaling evolution (see the 
specific example in Eq. (p?7)l below). 

Interestingly, the evolution of the scaling parameters 
6i, and therefore the nontrivial line element H15|) . can be 
generated already in a non-expanding cloud with time- 
dependent interaction g{t). A similar experiment has 
been suggested in jsj], where time dependent interac- 
tions were used to simulate FRW cosmologies and quan- 
tum quasiparticle production. The difference to our ap- 
proach is due to the fact that the authors of [s^l con- 
sider a trap with very steep walls (effectively a hard- 
walled container), so that the density of the cloud does 
not change, and the superfluid velocity vanishes every- 
where at all times. In our setup, we are able to induce 
cosmological quasiparticle production in a harmonically 
trapped gas, by changing simultaneously the harmonic 
trapping and the interaction. The simplest case we can 
consider is to leave all 6^ = 1, like in a stationary Bose 
condensate. We then create "cosmological" quasiparti- 
cles just by changing g (using Feshbach resonances, cf.. 
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e.g., Refs 
cies uji = 
(PI and im 



[13 ) , and accordingly change the trap frequen- 
ts (in the isotropic case). Following Eqs. (jS)), 



we then have the simple relations 



^^0 5(0) 



F{ty 



(32) 



The metric associated with such a thermal quasiparticle 
universe created by "shaking the trap" and simultane- 
ously changing the interaction appropriately reads, from 
Eq. ^ 



.9(0) 



(33) 



Now, defining the scale factor of the BEC quasiparticle 
universe by the relation 



seal 



5(<)/5(0)' 

we have, up to the (irrelevant) factor l/g{0), 
ds^ = -a^^^idr^ -t- dx^^. 



(34) 



(35) 



This is the form of the metric employed in '34'|, where it 
was used to calculate cosmological quas iparticle produc- 
tion, inspired by a model of Parker |^ . Note that here 
a nontrivial scale factor is induced without expanding 
the cloud. We see from relations (|32|l and H34|l that the 
scale factor Ogcai is in our harmonically trapped case sim- 
ply proportional to the ratio of initial and instantaneous 
trapping frequencies, Cscai = Vcuoo/ijj{t). 

In Ref. |3j| , a specific choice of the scale function 
o.sca\{Ts) was taken for the calculation of the quasipar- 
ticle creation process. 



-^scal 



^scal,i + "scal.f 



cal.f 



■ tanh 



where Oscai.i 
respectively. 



(36) 

and ascai.f are initial and final scale factors, 
In the adiabatic approximation, one ob- 



latic ap 

tains a thermal spectrum j34l l35l |. with a temperature 
governed by the inverse laboratory time scale to ^ ^so 
on which trapping frequencies, interaction and thus the 
scale factors change: 



Tn 



1 



scal.i 



■ '^scal.f 



47rto a^, 



(37) 



This temperature is, according to Eq. H30() . determined 
by the singular points of the tanh function in Eq. H36II . 
The fact that the spectrum is thermal is obtained in J|4J 
for a specific example with a certain form of the time 
dependent interaction. We emphasize here that the ther- 
mal spectrum is a generic feature of adiabatic evolution 
in harmonically trapped superfluid gases with temporally 
varying trapping potential and interactions. 



V. DETECTION BY MEASURING 
DENSITY-DENSITY CORRELATIONS 

Although the solutions of the hydrodynamic equations 
are unique, their interpretation in terms of the number of 
phonons in a given mode is subject to all the conceptual 
difficulties encountered by the definition of particle states 
in curved space-times 0. In [2& |. we have shown that 
simply by choosing a specific realization of a quasiparticle 
(phonon) detector one can observe thermal quantum "ra- 
diation" from a de Sitter horizon (the Gibbons-Hawking 
effect 32] ) as a purely choice-of-observer related phe- 
nomenon, without energy transport or dissipation taking 
place inside the liquid. Below, we confine ourselves to 
the standard (conventional) laboratory means of parti- 
cle detection (a CCD camera detecting individual atoms, 
rather than phonons), and concentrate on uniquely de- 
fined laboratory frame observables, such as the lab frame 
density-density correlations discussed in what follows. 

The density fluctuation operator is given by 



5p - E 



1 



d 



2Venk dTs 



so that the lab- frame density-density correlator G{xi2) — 
{6p{xi)Sp{x2)) /V"^ , averaged over the initial state is, in 
the isotropic case. 



F 



2Vk V2 ' 
[l + 2|/?„ 



2Re 



(39) 



-2ie„ / d; 



Here the normalization condition H28|l is used, and ri2 = 
|a3i2| = \xi - X2\ < |a;i|, \x2\. 

The TP condition H26|l ensures that the cross-term 
proportional to af3* averages to zero at large times t. 
The term with 1 in the square brackets describes the 
evolution of the vacuum fluctuations and the summa- 
tion over n is cut off at the Planckian energy scale: 
max[e„] ^ Epi — pg{0). Accordingly, the correspond- 
ing correlation function decays at Planckian distances 
ri2/6 ^ c/Ep\ and is very short-range. Subtracting the 
vacuum contribution, we obtain the following expression 
for the regularized correlator 



Grog(a;i2) 



We note that in QFT the regularization procedure does 
not follow in a unique manner from the field theory it- 
self, and can be applied using different assumptions about 
the high-energy behaviour of the excitations created from 
the fundamental "ether." Here, the spectrum (and ori- 
gin) of the TP excitations is well-known, and hence the 
above regularization of two-point correlation functions 
can always be strictly justified. This regularization pro- 
cedure of course is not limited to density-density correla- 
tors only. A similar technique can be used, for example, 
to find a regularized energy-momentum tensor. 
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To be more specific, consider a large Bose-Einstein- 
condensed gas cloud in the Thomas- Fermi hmit. Then, 
close to the center of the gas we can use WKB (plane 
wave) functions (/>„, with energies = cfc, and the regu- 
larized Green function is given by 

G..eg(ri2) - (^) , (41) 

where the function 

1 f°° sin(fcr/6) , , , 

The function Q reaches its maximum for ri2 = 0, so that 
the signal to noise ratio is maximally 

The above discussion shows that even from a small 
"noise" signal, one can extract the relevant features of 
the quantum state of the gas cloud (for a more detailed 
discussion cf. [13 )■ In order to be measurable, the 
quantity (|43|) has to be of the order of a few percent. 
This can in principle be achieved by using initially dense 
clouds with strong interparticle interactions. Finally, we 
mention that a similar, i.e., velocity- velocity instead of 
density-density "noise" correlation function has already 
been measured in the experiments of j37j| . 

VI. CONCLUSION 

In the present investigation, we have derived the gen- 
eral scaling equations for harmonically trapped super- 
fluid Bose and Fermi gases, and related these, in partic- 
ular, to quasiparticle metric tensors of the de Sitter and 
Friedmann-Robertson- Walker type, familiar from a cos- 
mological context. The quasiparticle creation in a har- 
monically trapped superfluid gas, by changing interaction 



and trapping simultaneously in an appropriate manner, 
can therefore be described in a general framework, and be 
interpreted to be analogous to the particle creation occur- 
ing during rapid expansion of the cosmos. In particular, 
it was found that for a readily experimentally available 
case, the harmonically trapped, dilute superfluid Bose 
gas, a FRW type metric can be induced if trapping fre- 
quency to and interaction coupling g are changed such 
that iv^{t) oc g{t), without expanding the gas. The cos- 
mological scale factor in this case is inversely proportional 
to the trap frequency, ascai(*) cx l/w(i). 

If the frequency mixing leading to quasiparticle cre- 
ation can be described in the WKB approximation, gen- 
erally a thermal distribution is found, where the temper- 
ature is determined by the singular points of the scaling 
factors given by Eq. (|10|) , in the complex plane of scaling 
time Tg. 

We finally stress that, in contrast to a typical cosmo- 
logical calculation, hydrodynamic fluctuations in a labo- 
ratory experiment always have a well-defined initial state 
in the lab frame, with time coordinate t. Therefore, am- 
biguities of the final quantum state as regards the de- 
pendence of its particle content on the initial conditions 
imposed on the "vacuum" can be ruled out: There exists 
the preferred lab frame vacuum, uniquely prescribing the 
initial particle content of the quantum field. 
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